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We present a concept of non-Gaussian measurement composed of a non-Gaussian ancillary state, linear optics 
and adaptive heterodyne measurement, and on the basis of this we also propose a simple scheme of implementing 
a quantum cubic gate on a traveling light beam. In analysis of the cubic gate in the Heisenberg representation, 
we find that nonlinearity of the gate is independent from nonclassicality; the nonlinearity is generated solely 
by a classical nonlinear adaptive control in a measurement-and-feedforward process while the nonclassicality is 
attached by the non-Gaussian ancilla that suppresses excess noise in the output. By exploiting the noise term as 
a figure of merit, we consider the optimum non-Gaussian ancilla that can be prepared within reach of current 
technologies and discuss performance of the gate. It is a crucial step towards experimental implementation of 
the quantum cubic gate. 

PACS numbers: 03.67.Lx, 42.50.Dv, 42.50.Ex, 42.65.-k 


I. INTRODUCTION 

Development and application of quantum physics crucially 
rely on progress in quantum operations with various physi¬ 
cal systems. For discrete-variable systems, a basic controlled- 
NOT nonlinear gate [1] has been already demonstrated with 
many systems [2-5] and the current problem is scalability of 
their implementations. On the other hand, for more complex 
continuous-variable (CV) systems [6], a full set of basic oper¬ 
ations has not been closed yet. It was proven that in order to 
synthesize an arbitrary unitary operation, it is enough to add 
a cubic nonlinear operation to the already existing Gaussian 
operations [7]. Any nonlinearity can be principally obtained 
from a chain of the Gaussian operations, the cubic nonlinear¬ 
ity and feedforward corrections [7, 8]. The cubic nonlinearity 
is therefore a bottleneck of CV quantum physics. 

Already a decade ago, Gottesman, Kitaev and Preskill 
(GKP) suggested a way how to implement a cubic nonlinear 
gate based on Gaussian operations, Gaussian measurement, 
quadratic feedforward correction and an ancillary cubic state 
produced by the cubic nonlinearity [9] . Various approaches to¬ 
wards the cubic gate have followed [10-13]. Particularly in the 
field of quantum optics, most of the components of the cubic 
gate have been experimentally demonstrated, mainly because 
of the high quality of generating squeezed states and efficient 
homodyne detection. The Gaussian operations have been al¬ 
ready mastered [14-16], utilizing a concept of measurement- 
induced operations [17]. Furthermore, they have been tested 
on non-Gaussian states of light [18], to prove their general ap¬ 
plicability. Recently, the quadratic electrooptical feedforward 
control has been demonstrated [19]. In addition, to indepen¬ 
dently obtain the cubic state, a finite dimensional approxima¬ 
tion of the cubic state has been suggested [20] and its perfor¬ 
mance in the GKP scheme has been discussed. The cubic state 
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has been experimentally generated as a superposition of pho¬ 
tons and verified [21]. Potentially, such a superposition state 
can be stored in and retrieved from recently-developed opti¬ 
cal quantum memories [22, 23]. In order to make resource 
nonclassical states compatible with the measurement-based 
scheme, real-time quadrature measurement of a single-photon 
state has been demonstrated [24]. 

A drawback of the original GKP idea is that requires to 
implement the quantum nondemolition gate, i.e., the CV 
controlled-NOT gate [17], and a squeezing feedforward that 
depends on the measurement result. While each of them has 
been already demonstrated [15, 19], the total implementation 
to build a unitary cubic operation demands three squeezed 
states as well as one non-Gaussian ancilla, and is probably not 
the simplest arrangement. In contrast, we here use adaptabil¬ 
ity of linear optical schemes and propose a better and simpler 
topology with linear optics and suitable ancillary states. 

Our approach is to tuck all the non-Gaussian aspects into 
the measurement process. The topology will be then similar to 
the simple one used for a measurement-induced squeezing gate 
[14, 17-19, 25]. Non-Gaussian operations can be realized by 
simply substituting a measurement of nonlinear combination 
of quadrature amplitudes for the Gaussian homodyne measure¬ 
ment [26, 27]. We construct such a measurement in a form of 
a generalized non-Gaussian measurement by combining ordi¬ 
nary Gaussian measurement tools with non-Gaussian ancillary 
states that can be prepared with photon detection. In fact, we 
can exploit arbitrary superpositions of photon-number states 
up to certain photon level within reach of current technologies 
[21,28,29]. 

In this paper, we first provide an idea of non-Gaussian mea¬ 
surement comprising a non-Gaussian ancillary state, linear op¬ 
tics and adaptive heterodyne measurement. Using the non- 
Gaussian measurement, we next propose a simple schematic of 
a quantum cubic gate based on the measurement-induced oper¬ 
ation scheme, whose resource states are only one squeezed vac¬ 
uum and one non-Gaussian state. While in previous work the 
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FIG. 1. (Color online) Minimal single-mode implementation of 
measurement-induced quantum operation. 


input-output relation of the cubic gate has been investigated in 
the Schrodinger picture, here we analyze the gate in the Heisen¬ 
berg picture to include imperfections in the scheme. We then 
find that nonlinearity of the gate is independent from nonclas- 
sicality. Specifically, the nonlinearity is generated solely by 
a classical nonlinear adaptive control in a measurement-and- 
feedforward process regardless of the non-Gaussian ancilla. 
On the other hand, the nonclassicality is attached by the an¬ 
cilla that compensates residual noise in the output. Finally, we 
discuss an overall performance of the cubic gate in such a topol¬ 
ogy and consider non-Gaussian ancillary superposition states 
up to certain photon level to investigate how well the unwanted 
noise can be suppressed in the gate. 


II. MINIMAL IMPLEMENTATION OF 

MEASUREMENT-INDUCED QUANTUM OPERATIONS 

Measurement-induced quantum operation scheme [17] de¬ 
composes various quadratic operations into linear optics, dis¬ 
placement operation, homodyne detection and offline squeezed 
light beams, which are readily available in actual optical ex¬ 
periments. One of the realizations of the scheme is the basic 
squeezing gate. Firstly we combine an input state |t/)) and an 
eigenstate \x = 0 ) of the position quadrature x at a beam split¬ 
ter whose transmittance is represented by y/T. We then mea¬ 
sure the momentum quadrature p of one of the optical modes 
and obtain a value y. Finally we apply displacement to the 
p quadrature of the remaining mode with the value Pdisp = 
y/(l — T)/Ty and obtain a squeezed output state. Ideally the 
output is a pure state S\^p) where S is an x-squeezing opera¬ 
tor defined as S^xS = yfTx and S^pS = p/y/T. In the case 
of implementing p squeezing, it is enough to replace the ancil¬ 
lary x eigenstate with the p eigenstate \p = 0 ) and exchange 
the roles of x and p quadratures. This type of operation has 
been successfully demonstrated in [14, 18], where the position 
eigenstate is replaced with the squeezed vacuum. 

On the basis of one-way CV cluster computation [26, 27], 
we can generalize the basic squeezing gate to minimal single¬ 
mode implementation of arbitrary-order quantum operations as 
shown in Fig. 1. The homodyne detector in the squeezing gate 
is now replaced with a detector that measures a general quadra¬ 
ture U^(x)pU n (x), where the unitary operator U n (x) is de¬ 


fined as nth-order phase gate U n (x) = exp(i 7 x n ) with a real 
parameter 7 . Hereafter we set h = 1 for simplicity. The mea¬ 
sured general quadrature is thus p + njx n ~ 1 . In the ideal case, 
the output state is expressed as SU n (y/ 1 — Tx)\i/j). This gate 
deterministically applies the phase gate to the input state with 
the additional constant squeezing that can be compensated by 
another squeezer. 

It is known that, arbitrary single-mode unitary can be de¬ 
composed into the set of gates U n (x) for n = 1, 2 ,3 for all 
7 G M, together with the 7t/2 phase shift [8, 27]. This also 
holds when we exploit the minimal implementation in Fig. 1 . 
Ui(x) is the trivial displacement operation, and 62 (x) has been 
experimentally demonstrated [19, 25]. The remained task is 
thus to realize a cubic gate Us(x). We now consider how to 
construct measurement of the nonlinear quadrature p + 37 X 2 
with affordable apparatuses, as explained in the following sec¬ 
tions. 


III. NON-GAUSSIAN MEASUREMENT BY GENERALIZED 
HETERODYNE DETECTION 

A. Projecting on pure states 

In quantum physics, measurements are represented by oper¬ 
ators. In the simplest case of von Neumann measurements, 
these operators are simply projectors on particular quantum 
states. In the case of the keystone measurement of CV quan¬ 
tum optics, the homodyne detection, each measurement result 
indicates the measured state was projected on an eigenstate of 
the measured quadrature operator. Analogously, the hetero¬ 
dyne detection, which can be modeled by a pair of homodyne 
detectors simultaneously measuring conjugate quadratures of 
a mode split by a balanced beam splitter [30], implements a 
projection onto a coherent state. Both of these kinds of mea¬ 
surements are Gaussian—the measured quadrature distribution 
is Gaussian if the measured state is Gaussian. 

One way to achieve a non-Gaussian measurement is to take 
advantage of non-Gaussian states in combination with the stan¬ 
dard heterodyne detection schemes. The basic idea of the mea¬ 
surement is best explained in the x representation. Consider we 
have a standard heterodyne detection configuration, where the 
idle port of the beam splitter is not injected by a vacuum but by 
a specifically prepared ancillary state |^a) = f^ a(x)\x) dx. 
For a particular pair of measurement results q and y , the pro¬ 
cedure implements projection onto a state 

D(V2q + iV2y)t\il>A)- ( 1 ) 

Here D(o ) = exp{z\/2a' Ira [a ] — iy/2pRe[ a] } stands for the 
displacement operator and T is the time-reversal antiunitary 
operator represented by T^xT = x and T^pT = —p. To 
derive the expression ( 1 ), we can start with the projection states 
of the pair of homodyne detectors 

(xi = q\(p 2 = y\. (2) 

If we take into account the unitary balanced beam splitter, the 
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projection state becomes 



q + x 2 

V2 


/ -q + x 2 

\ V2 


e ~iyx 2 ' 


( 3 ) 


During the measurement, this state will be jointly projected 
onto the measured and the ancillary state. The measured state 
is unknown, but we can already apply the ancilla in the second 
mode. This reduces the state to 


of delta functions S(x i — q)S(p 2 — y), which we then propa¬ 
gate through the beam splitter and apply to the ancillary state, 
resulting in 

Wu(xi,pi\q,y) = 

J dx 2 dp 2 w A (x2,P2)s ( i^j d ~yj • 

( 8 ) 



where the subscript was dropped because it was no longer 
needed. After a straightforward substitution we can express 
the projection state as 

J dx^ A {x)e l ^ xy \x + V2q). (5) 

Since the time-reversal operator corresponds to complex con¬ 
jugate in the x representation, the expression (5) is the same as 
Eq. (1). For q = y = 0, we obtain simple projection onto the 
given ancillary state f (x)\x) dx. We can see that if the an¬ 
cillary mode is in the vacuum or a coherent state, the measure¬ 
ment remains the simple heterodyne detection, as is expected. 
However, if the ancilla is non-Gaussian, we obtain a truly non- 
Gaussian measurement. 

B. Projecting on impure states 

In a realistic scenario, the ancillary state will be generally 
not pure. To take this into account, it is best to abandon the x 
representation and employ the formalism of Wigner functions. 
The basic premise, however, remains. The measurement still 
implements projection onto a specific state, only this time the 
state will be represented by a Wigner function. Specifically, 
for a signal two mode state represented by a Wigner function 
Ws (xq , po , x\ , pi ), the outcome of a measurement performed 
on mode 1 yielding a pair of values q and y results in the Wigner 
function 


C. Arbitrary Gaussian operations within the measurement 

One may desire to apply Gaussian operation to the non- 
Gaussian ancilla because some Gaussian operations (such as 
squeezing) enhance certain features of the state. Here we show 
that, instead of projecting on a raw non-Gaussian state, we can 
alter the measurement so it projects on a non-Gaussian state 
altered by an arbitrary Gaussian operation. This can be enor¬ 
mously useful because we do not need to implement an addi¬ 
tional Hamiltonian that often makes the state impure in actual 
experiments. Note that we are disregarding displacement be¬ 
cause that can be achieved simply by displacing the measure¬ 
ment results. For a pair of quadrature variables x and p, an ar¬ 
bitrary Gaussian operation is represented by a real two-by-two 
symplectic matrix S whose elements satisfy sns 22 — <^ 12^21 = 
1. If we consider that phase shift can be implemented “for 
free,” the arbitrary Gaussian unitary transformation reduces to 

x' = z\x, p f —p -\- z 2 x , (9) 

zi 

where z\ and z 2 are arbitrary real parameters. To achieve 
this transformation, we must modify the measurement setup 
in two ways. First, the balanced beam splitter in Eq. (3) will 
be removed and replaced by a beam splitter with transmittance 
T and reflectance R = 1 — T. Second, instead of measur¬ 
ing quadrature p 2 we measure quadrature p 2 (0) = p 2 cos 0 + 
x 2 sin 6. The projection functions of the measurements itself 
in Eq. (8) are then 

S(xi — q)S(p 2 cos 6 + x 2 sin# — y). (10) 


W OVLt {xo,po\q,y) 

f ( 6 ) 

oc / dx 1 dp 1 Ws{xo,Po,x 1 ,p 1 )W M (xi,Pi\q,y), 

where the function Wm (xi,Pi\q,y) represents the projector on 
the particular state. In our scenario, in which the pair of homo¬ 
dyne detectors are supplied with an ancillary state correspond¬ 
ing to a Wigner function Wa(x,p), the projector function can 
be found as 

Wu(x,p\q,y)—2 Wa(x - V2q,-p+V2y). (7) 

We can see that this form agrees with Eq. (1) if we realize that 
the time-reversal operator T transforms the Wigner function 
variables as (x,p) (x, —p). The relation (7) can be derived 

in the same way as relation (1). We start with the homodyne 
measurement projector functions, here represented by the pair 


Using the same steps we used to arrive at Eq. (7) we can now 
obtain the generalized projection function 


W M (x,p\q,y) = 


Wa | J I*- 


|/RTcos6>| ' 
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q 

Vr’ 


^p- 


tan 0 q tan 9 

^= x + 7= + 

y/RT VR 
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VT cos 6 ^ 
( 11 ) 


We can immediately see that after the time-reversal opera¬ 
tions, we have z\ = y/T/R and z 2 = tan 0/y/RT and these 
two parameters can attain arbitrary real values. As a conse¬ 
quence, after addition of a phase shift the function ( 11 ) imple¬ 
ments projection onto the ancillary state altered by an arbitrary 
Gaussian operation. 

It is worth pointing out that the two homodyne measure¬ 
ments need not to be independent. One of the measurements 
can have parameters changing based on the results of the other 
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one, thus creating a sort of adaptive measurement scheme. For 
example, the measurement phase 0 can depend on the measure¬ 
ment result q. This can be used to induce a nonlinear behavior, 
as we will see in Sec. IV B. 


IV. IMPLEMENTATION OF A CUBIC GATE 
A. With nonadaptive non-Gaussian measurement 

In this section we apply the non-Gaussian measurement to a 
particular task—the implementation of a nonlinear cubic gate 
U = e nx to an arbitrary quantum state. In terms of quadra¬ 
ture operators, the gate performs transformation 

x r = x, p'=p + 3'yx 2 . (12) 


Before proceeding to a scheme with the adaptive heterodyne 
measurement, we firstly consider implementation with non¬ 
adaptive measurement expressed by Eq. (1). 

The basic principle of the operation can be quickly explained 
in the x representation. The unknown input state |^) is mixed 
with a squeezed state on a balanced beam splitter. If we for ease 
of explanation consider the infinite squeezing, the resulting two 
mode state can be expressed as 


/ 


dx'ip(x) 




(13) 


After applying non-Gaussian measurement (4) on one of the 
modes, we obtain the projected state in the form 


Non-Gaussian measurement 



FIG. 2. (Color online) Schematic of a cubic gate. BBS, balanced 
beam splitter; HOM, homodyne measurement; LO, local oscillator; 
PS, phase shift; NL, nonlinear classical calculation. While all the 
optics are linear, the classical circuit involves nonlinear calculations 
that makes the feedforward nonlinear. The nonlinear classical circuits 
have been already devised in the experiment of dynamic squeezing 
[19]. 



(14) result. This allows the whole procedure to operate determinis¬ 
tically. 


where q and y are again the homodyne measurement results, 
and ^a{x) is the wave function of the ancillary state. For im¬ 
plementing the cubic gate, ipA (x) has to be cubically dependent 
on x , and is ideally in a state 

IV»a) = J dx ex.p(i'yx 3 )\x) (15) 

and the whole operation would lead to 


B. With adaptive non-Gaussian measurement 

In Eq. (16) we need quadratic feedforward in the form of ad¬ 
justable squeezing. Thus the topology here is not as simple as 
the minimal implementation depicted in Fig. 1 . To realize mea¬ 
surement of the nonlinear quadrature p + 3^/x 2 , we exploit the 
adaptive non-Gaussian heterodyne measurement. According 
to the results in Sec. IIIC, by altering the phase of the second 
measurement, we can project onto a transformed ancillary state 


exp(i7x 3 ) J dx!p(x) 


q 2 - V2y)x] 



V^q + if — a/2 q tan 0 ) 

\ cos 0 J 

X \ 

(16) 

D 

v^r 



L \ / J 


fe ^ 2tan0 |^ A ) 


(17) 


This is almost exactly the desired output state. The only dif¬ 
ference is a constant squeezing and two unitary operations de¬ 
pending on the measured values. The constant squeezing can 
be fully compensated either before or after the operation and 
the measurement dependent unitary operations can be removed 
by a proper feedforward. This is exactly the same principle 
as employed by the CV teleportation and CV measurement- 
induced operations. While each particular measurement result 
projects on a different quantum state, these states belong to the 
same family and the proper operation can smear the differences 
and produce a quantum state independent of the measurement 


Again, q and y are the measured values in the heterodyne de¬ 
tection, and IV’a) is the cubic state (15). We then substitute 
3x^7 q for tan 0. After simple algebras we find the projection 
state 


exp(— ijx 3 ) 



(18) 


which means an eigenstate of the nonlinear quadrature p+3'yx 2 
with the eigenvalue y/2yj cos 0. This scheme can be illustrated 
as Fig. 2. Here the quadrature basis of the second homodyne 
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detection is determined by the result of the first homodyne de¬ 
tection. As a result, the heterodyne detection and the classical 
calculation compose a module of non-Gaussian measurement, 
and the feedforward is now the simple displacement operation. 

To explicitly show how this scheme works, it is instructive to 
employ the Heisenberg representation, which would have the 
added benefit of incorporating the imperfections arising from 
the realistic experimental implementation, e.g. finite squeez¬ 
ing. May the unknown signal mode be labeled by ‘O’ and de¬ 
scribed by quadrature operators xo and po. After combining 
the initial state in mode ‘ 0 ’ with the squeezed state in mode 4 1 ’ 
and with the non-Gaussian ancilla in mode ‘2’, the respective 
quadrature operators read 




Xi = -Xo + -Xl - -J=X 2 , 

p[ = \po + \h - ^ 2 , 


V V 1 . 

x 2 = 2 X ° + 2 Xl + 

& = \*> + \h + 


(19a) 

(19b) 

( 20 a) 

( 20 b) 

( 21 a) 

( 21 b) 


In the next step we measure the x quadrature of mode V 
and obtain value q. We can now use the value to adjust the 
measured phase of the second homodyne detector. In ef¬ 
fect we end up measuring the value y of quadrature operator 
^2 sin 0 -f p 2 cos 0, where 0 = arctan( 3 \/ 27 </). Note that, 
since 0 nonlinearly depends on q, which carries information of 
x\ quadrature, we can interpret this type of measurement as 
the origin of nonlinearity of the gate. The quadrature opera¬ 
tors of the output mode can be now expressed in terms of the 
measured values as 


Po = V2po+P2 + y[(io + Xl ) 2 - 2x 2 ] 


( 22 a) 


ypv_ 

cos 0 


( 22 b) 


The last term of the p quadrature, which is the only term explic¬ 
itly depending on the measured values q and y, can be removed 
by a suitable displacement and we are then left with the final 
form of the operators: 


p" = x/2 (po+2^r 2 ) 

+ (p 2 - 37 X 2 ) + 37 (x 0 xi + if 2 


(23a) 


(23b) 


Both of the first terms in Eq. (23) represent the ideal cubic op¬ 
eration, i.e. combination of the cubic gate e l i(x 0 /V2) an( j the 
constant squeezing mentioned in Sec. II. Those terms does not 
depend on the quadratures of the other ancillary states. Differ¬ 
ently from the output (16) in Sec. IV A, in the Heisenberg rep¬ 
resentation we can say that the cubic nonlinearity comes from 
the adaptive non-Gaussian measurement and feedforward pro¬ 
cess regardless of the ancillary states. 

Naturally, the ancillary states are still required to complete 
the operation since the outputs have residual terms. It is 
straightforward to find the ideal ancillary state in the mode 4 1 ’ 
the quadrature eigenstate \x = 0)i because the state affects 
only on the last terms of Eq. (23) and they vanish when x\ 0. 
In experimental implementation, we approach the ideal state by 
using squeezed vacuum states. 

On the other hand, the middle term of Eq. (23b), Pnlq = 
P 2 — 37 x 2 , depends solely on the ancilla in the mode ‘2’. This 
term vanishes when the ancilla is the ideal cubic state (15). 
This state is best approached by considering physical states 
that squeeze the nonlinear quadrature Pnlq, as discussed in 
the next section. 


V. OPTIMUM ANCILLARY STATE 

To find suitable states in the mode ‘2’, we can use the ex¬ 
pectation value and the variance of the nonlinear quadrature 
Pnlq as figures of merit, both of which should be approach¬ 
ing zero. Here we consider preparing the ancillary state that 
can be generated within reach of current technologies. On one 
hand, arbitrary superpositions of photon-number states up to 
three photon level |Vw= 3 ) can be prepared [21, 29], and the 
photon-number limit can be in principle incremented. On the 
other hand, we can perform universal Gaussian operation Uq 
onto any input state [14-16]. Then the ancilla best suited for 
our purposes can be found in a form UgI'Pn) by optimizing 
over all superposition states up to 7V-photon level \P>n) and all 
Gaussian operations Uq that can be applied on the state after¬ 
wards. In this way, we are using the expensive non-Gaussian 
resources only for the key non-Gaussian features of the state 
[3 i] . 

Our goal is to find a state UgI'Pn) that minimizes the expec¬ 
tation value (pnlq) and the variance V(pnlq) = ((pnlq - 
(Pnlq)) 2 )- The operator is symmetric with respect to space 
inversion, x 2 —)> — £ 2 , and has a linear term of P 2 . Accord¬ 
ingly the relevant Gaussian operations are the p displacement 
represented by p 2 P 2 and the x squeezing represented 
by X 2 X 2 /X and P 2 Xp 2 , where po and A are arbitrary 
real parameters. Thus the nonlinear quadrature after suitable 
Gaussian operations is represented as 


uIpnlqUg = 7 1/3 


A'P2 —3(^ 


2 


+ P 0 , 


(24) 


where A' = A/ 7 1 / 3 . From this point of view, we can see that 
the expectation value (pnlq) vanishes when we apply suitable 
displacement po- On the other hand, the variance V (pnlq) can 
be minimized by optimizing the state \P>n) and the parameter 
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FIG. 3. (Color online) Variances of the nonlinear quadrature with the 
optimized photon-number-state superpositions up to N photons, nor¬ 
malized by the Gaussian limit The parameter X' is optimized 

over to find the minimum of the variance. 



FIG. 4. (Color online) Absolute values of the coefficients of the opti¬ 
mal finite approximation of ancillary states for various upper bounds 
of photon number N. Note that the coefficients of even-number pho¬ 
tons are real and odd-number photons imaginary, due to the symmetry 
of the nonlinear quadrature p — 3jx 2 with respect to x —> —x. 


A'. Furthermore, since A' can be any real number, we can say 
that the optimum state does not depend on 7 . We therefore 
use the variance of X'p 2 — 3(£ 2 /A') 2 as the actual figure of 
merit to derive the optimum state |Vl/v) and the corresponding 
parameter A'. 

Let V£ pt the minimum value of the variance V (pnlq) with 
the optimum state | 7 / 7 ^) anc * the optimal parameter A /opt . 
Note that Vo° Pt represents the Gaussian limit—the minimum 
variances when the state is optimized over all Gaussian states. 
Then the relative noise V^ pt /VJ)° pt , as shown in Fig. 3, repre¬ 
sents the ratio of the minimum noise to the Gaussian limit, and 
is independent from 7 . We can see that the variance decreases 
approaching zero with N, and that even a state obtained as a 
superposition of zero and one photon gives a substantial bene¬ 
fit over the Gaussian limit. To present the optimized states, we 
represent the optimal approximate state up to 7V-photon level 
by | = En=o c n Pt | n ) an d pl°t absolute values of the co- 



X X 


FIG. 5. (Color online) Wigner functions of the optimal ancillary 
states, (a) The ideal cubic state for 7 = 0.1 (normalized over the 
displayed area), (b) N = 1, (c) N = 3, (d) N = 5, (e) N = 9. Note 
that the approximate states have offsets in the p direction, which can 
be compensated by p displacement. 


efficients in Fig. 4. In the case of optimizing the superposition 
state up to three photons, the optimal approximate state looks 
as 


1 ^ 3 } OC 0.17|0) - 0.56z|l) - 0.73|2) + 0.35z|3). (25) 


The state is different from the cubic state from [20] because 
of the different derivation of the states. In [20], the state was 
determined as if it was produced by applying the cubic gate to 
the vacuum without considering optimization over squeezing 
and displacement. On the other hand, the present state (25) is 
derived so that its overall suitability as the ancilla is maximized 
with suitable Gaussian operations. In either case, the state can 
be prepared by the same experimental method [21, 29]. 

The cubic nature of the states is also nicely visible from their 
Wigner functions as depicted in Fig. 5. For comparison, we 
check the Wigner function of the ideal cubic state [11] 


W (x,p) = 2nJ\f 


4 

1/3 . ( 

" 4 " 

37 

A1 

.37. 


1 1/3 


[ 37 a : 2 - p) , (26) 


where Ai(x) is the Airy function and AT a normalization factor. 
As can be seen from Fig. 5 (a), the Wigner function is sym¬ 
metric with respect to the p axis and has a oscillating parabolic 
shape. These characteristics also appear in the approximate 
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cubic states shown in Fig. 5 (b)-(e). As the upper limit of pho¬ 
ton number becomes larger, the number of fringes along the p 
direction increases approaching the ideal one. Those Wigner 
functions of the approximate states can be considered to show 
core non-Gaussianity that then spreads out on the phase space 
by the following optimized squeezing. 

So far we have not consider how to implement the optimized 
squeezing onto the core non-Gaussian state. Actually, instead 
of adding another squeezing gate, the squeezing operation can 
be embedded into the adaptive non-Gaussian measurement by 
using the results in Sec. IIIC. We discuss the details on it in 
Appendix A. 


VI. CONCLUSION 

We have introduced the concept of an adaptive non-Gaussian 
measurement—a CV measurement with a set of possible val¬ 
ues, each of which is associated with a projection onto a non- 
Gaussian state. The measurement is realized by a pair of homo¬ 
dyne detectors and a supply of suitable non-Gaussian ancillary 
states. One particular advantage of this measurement is that an 
arbitrary Gaussian operation can be implemented on the soon 
to be measured quantum system simply by tools of passive lin¬ 
ear optics. In addition, some non-Gaussian operations can be 
implemented in a same way by making some of the measure¬ 
ment parameters dependent on already measured values. 

To demonstrate this design feature, we have proposed a new 
method of realizing the cubic gate [20]. The current proposal 
does not require active operations to be performed on the trans¬ 
formed quantum system, all of them being part of the non- 
Gaussian measurement, which significantly improves the fea¬ 
sibility of the setup. Specifically in the Heisenberg represen¬ 
tation, it turns out that nonlinearity of the gate is created clas¬ 
sically while the nonclassicality is given by the non-Gaussian 
ancilla in terms of reducing residual noise. By exploiting the 
noise term as a figure of merit, we have found a new class of 
ancillary states that promise better performance than the states 
of [21]. The final implementation of the complete cubic gate 
can be therefore expected soon. 


We explain here an input-output relationship of the cubic 
gate with unbalanced beam splitters. Transmittance and re¬ 
flectance of the first beam splitter are represented as T\ and 
R\ = 1 — T\ , respectively. T 2 and R 2 are also defined in the 
same way for the second beam splitter. After the beam splitter 
transformations, the quadratures of the output modes are 


Xq = \ff\x0 - \/-R/ 2 x\, 

(Ala) 

0 " 

II 

0 

1 

to 

(Alb) 

x\ = \JR\T 2 xq + \/T{T 2 x\ - \J-R 2 x 2 , 

(A2a) 

Pi = \J R 1 T 2 P 0 + \/T\T 2 pi — \^R 2 p 2 ^ 

(A2b) 

x' 2 — \/ RiR 2 xq + \J T\R 2 x\ + 's/t^x 2 , 

(A3 a) 

P 2 = V R 1 R 2 P 0 + \/TiR 2 pi + \fT 2 p 2 . 

(A3b) 


After measuring the x quadrature of mode 1' and obtain value 
q , we set the phase factor 


0 = arctan 


6 T 27 \ 

VRiV- 


(A4) 


Then we measure the quadrature x 2 sin 6 + p 2 cos 0 and obtain 
value y. The p quadrature of the unmeasured mode O' can be 
expressed with the measured values q and y as 


Po = 


VR 


+ 


Po - 




RlT 2 

RR 2 


VTJhcose 

R±(R \ 3/2 


P 2 - 671 


R \R 2 J 


( 6i?iT 2 7 „ 6\/R[T 2 ^ 


\VR 


? 3/2 


%0 + 


VT\r 2 ± t 2 


r: 


, 3/2 


-xi q. 


(A5) 


We apply p displacement to this quadrature with value 
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Appendix A: Cubic gate with unbalanced adaptive 
non-Gaussian measurement 

By replacing BBS2 in Fig. 2 with an unbalanced beam split¬ 
ter, we have another degree of freedom to effectively apply ar¬ 
bitrary squeezing operation onto the ancillary non-Gaussian 
state, as shown in Sec. IIIC. Thus the Gaussian optimiza¬ 
tion discussed in Sec. V can be embedded in the cubic-gate 
schematic. 


Pdisp — 


VRl 


;2/ + 


^VRiT 2 (T 2 — R 2 ) 2 


VR R ‘2 cos 6" s/r\R 

and obtain the output quadratures 


, 3/2 

'2 


A" 

Xq — 


Po = 




VR ( *0 - 






RiT 2 \ 

r 2 ) 


3/2 


RlT 2 

r 2 


(P2 - 37*2) 

3/2 


+ 


67S,V^(g) (io ii + | 


(A6) 


(A7a) 


-x\ 


(A7b) 
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We can see that the outputs are equal to Eq. (23) if we set T\ = 
Ri = T 2 = R 2 = 1 / 2 . Note that, if we use unbalanced beam 
splitters, the displacement has a quadratic term as shown in 
Eq. (A 6 ). 

To explicitly see how the transmittances of the beam split¬ 
ters affect on the quadratures of the ancillary non-Gaussian 
state, we scale the strength of cubic nonlinearity 7 to 
(R 2 /R 1 T 2 ) 3 / 2 The output p quadrature (A7b) is then ex¬ 
pressed as 


i/(l^),A') = (V’|[y(A , )-(y(A , ))] 2 l^), (Bib) 

where (y( A')) is the expectation value (t/|?/(A / )|' 0 ). 

To make this problem digestible, we alternatively consider 
another minimization problem. Let d be a real number. We 
then replace the expectation value in Eq. (B lb) with d and set 
a new evaluation function 

Z(\iP),\',d) = (iPMX)-d} 2 M. (B2) 


Po = 


1 


vA 
+ 


C PO + 3 7 *o) 


RiT 2 

Ri 


P2 ~ 37 


R2 

RlT 2 


X2 


Next, we introduce another evaluation function W(d) defined 
as minimum of Z(|^), A',d) with respect to | ip) G Hn and 
A' G M. This can be expressed as 


(A 8 ) W(d) = min Z(\ip), A', d). (B3) 

v 7 Wen N 

x' gik. 



The second term represents the nonlinear noise determined by 
the non-Gaussian measurement. We can see that the ancilla 
is effectively squeezed by the squeezing factor yjR\T 2 /R 2 , 
which can be fully controlled by choosing transmittance of the 
second beam splitter. While universal squeezing operation in 
actual experiments [14, 18, 19] adds nonnegligible noise to the 
input state because of finite squeezing in its resource state, the 
effective squeezing in the heterodyne measurement does not 
require additional resource states, which helps preparation of 
the approximate cubic state with high purity. 


Suppose W (d) is minimum when d = d*. In addition, suppose 
Zd'ip), A', d*) is minimum when \^) = |^*) and A' = A'*. 
Then we can say that the set (1^*), A'*) is the true optimum 
set that minimizes V(\i/j),X'). This is verified as follows. Let 
(y(X'))* be the expectation value (V’*|f/(A / )|'0*). Then 




(B4) 


and consequently ((#(A'*))* — d*) 2 < 0 , which means d* = 
(#(A'*))*. Therefore for any \^) G Hn , any A' G M and 
the corresponding expectation value (y( A')) = {^\y{\')\^), 
it holds that 


Appendix B: Numerical method of approximating 
photon-number superposition to the cubic state 

In Sec. V, we considered the variance I^(pnlq) as a figure 
of merit to approximate the cubic state with photon-number- 
superposition states up to certain photon level and squeezing. 
Intuitively, the approximation can be done by numerically op¬ 
timizing all of the coefficients of a superposition state and the 
squeezing level, but it often leads to locally optimum solutions, 
especially when increasing the upper limit of photon numbers. 
Here we reduce the problem into optimization with two vari¬ 
ables, regardless of the size of the Hilbert space. With each set 
of the two variables, an optimized superposition state can be 
derived as an eigenstate of the minimum eigenvalue of a cer¬ 
tain positive-semidefinite operator. By numerically creating a 
minimum-search map with the two variables, we can make sure 
that the solution is almost certainly the true optimum one. The 
method is a variation of the classical variance-minimization 
problem [32]. 

Suppose Hn is a (N + 1)-dimensional Hilbert space up to 
7V-photon level, and |-0) is a state in Hn- Our purpose is to 
find a set of the optimum state \^) and the optimum parame¬ 
ter X' that minimizes the variance of the nonlinear quadrature 
y(X') = X'p — 3(x/X') 2 . This problem can be written as 

min VM),X'), (Bla) 

A G®. 


y(|V>*),A'*) = W(d *) 

< wm a '))) < (mm = a (B5) 

which means A'*) is minimum. As a result, the prob¬ 

lem can be solved by searching for a state that minimizes 
Zd'ip), A', d) with every A' and d. 

The point is that Z( , A', d) is a quadratic form, and there¬ 
fore each optimum state is determined as an eigenstate of the 
minimum eigenvalue of [y(X') — d] 2 represented by the limited 
Hilbert space. In the case that we look for the optimum state 
up to 7V-photon level, the matrix representation of [y(X') — d] 2 
reads 

N 

Y(X',d) = Y, Ymn(X,d)\m)(n\, (B 6 ) 

m,n=0 

Y mn (X',d) = (m\[y(X')-d} 2 \n), (B7) 

and the optimum state in terms of (A', d) is found as the eigen¬ 
state of the minimum eigenvalue of the matrix Y(X',d ), which 
can be deterministically obtained by numerical calculation. 
This implies that the problem is now broken down into a two- 
variable optimization problem. We can create a minimum- 
search map Z(\xp), A', d) with respect to A' and d , 

which makes it easy to look for the true optimum solution. 

Figure 6 shows some examples of the map used to derive the 
optimized superposition states in Fig. 5. We can see that the 
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-10 -5 0 -10 -5 0 

d d 



-10 -5 0 -10 -5 0 

d d 


number of local minima increases as the upper limit of photon 
numbers becomes larger. By choosing suitable ranges and res¬ 
olutions of (A', d), we almost certainly find the true minimum 
and consequently the true optimized state. 


FIG. 6 . (Color online) Minimum-search map. (a) N = 1, (b) N — 3, 
(c) N = 5, (d) N = 9. The values are normalized by shot noise level 
and shown in dB scale. 
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